On the linear response and scattering of an interacting molecule-metal system 
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A many-body Green's function approach to the microscopic theory of plasmon-enhanced spec- 
troscopy is presented within the context of localized surface-plasmon resonance spectroscopy and 
applied to investigate the coupling between quantum-molecular and classical-plasmonic resonances 
in monolayer-coated silver nanoparticles. Electronic propagators or Green's functions, account- 
ing for the repeated polarization interaction between a single molecule and its image in a nearby 
nanoscale metal, are explicitly computed and used to construct the linear-response properties of the 
combined molecule-metal system to an external electromagnetic perturbation. Shifting and finite 
lifetime of states appear rigorously and automatically within our approach and reveal an intricate 
coupling between molecule and metal not fully described by previous theories. Self-consistent incor- 
poration of this quantum-molecular response into the continuum-electromagnetic scattering of the 
molecule-metal target is exploited to compute the localized surface-plasmon resonance wavelength 
shift with respect to the bare metal from first principles. 



Ever increasing experimental interest in a variety of 
plasmon-enhanced molecular spectroscopies has provided 
impetus for the development of a corresponding as- 
sortment of theoretical descriptions of these phenomena 
[H, [3, H, 0, H, Si . Linear molecular spectroscopies such 
as Raman and fluorescence have found their plasmon- 
enhanced analogs in surface-enhanced Raman scatter- 
ing (SERS) and surface-enhanced fluorescence; both are 
now routinely realized in the extreme limit of single- 
molecule detection; see, e.g., Refs. i, i, And 
plasmon-enhanced versions of n-wave mixing and hyper- 
Raman scattering UlllJl, as well as other nonlinear spec- 
troscopies, are rapidly being explored as ultrasensitive 
probes of molecular structure complementary to those 
linear. Conversely, spectroscopy of the plasmon itself, 
either localized to metal particles or propagating across 
metal surfaces 13|, and embedded within a potentially 

0, 




resonant molecular environment [14l . Il5l | forms the basis 
for yet another promising direction of related research. 
Outside of the laboratory, the challenge remains to de- 
velop theories that are rich enough to describe the basic 
physics relevant to each, yet are simultaneously practi- 
cal enough to implement numerically for real systems of 
current experimental interest. 

Here we present the first numerical implementation of 
a new and general many-body Green's function formal- 
ism that is capable of describing a variety of plasmon- 
enhanced linear spectroscopies and apply it to model re- 
cent localized surface-plasmon resonance (LSPR) spec- 
troscopy and molecular sensing experiments. LSPR spec- 
troscopy is highly sensitive to the small refractive-index 
changes of the environment surrounding nanoscale metal 
particles and has been reported to detect the presence of 
as few as tens to hundreds of nearby molecules 
This sensitivity arises from the resonant coupling of inci- 
dent light to the collective oscillation of conduction elec- 
trons in the metal modified by the presence of the sur- 
rounding medium, which may have its own distinct elec- 



FIG. 1: Polarization effects between a quantum-mechanical 
molecule and a classical nanoscale metal are built into the in- 
teracting one-body molecular Green's function Q through the 
electron-plasmon potential Vp. Truncation of the associated 
self energy E* ~ E*j[Q:p(ep)] at second order in perturba- 
tion theory, followed by a partial resummation of the Dyson 
series for Q, allows the molecule to repeatedly polarize the 
metal's conduction electrons to all orders in perturbation the- 
ory. Underlying this interaction is the metallic polarizability 
Q.p{ui), which is connected to the macroscopic permittivity 
ep(k = 0,cij) of the metal through a particular model of di- 
electric response. 



tronic resonances. For metal spheres small in comparison 
to the wavelength of light, this resonance condition would 
be met at those frequencies w» in which 



Rc{ep(cj)/eA/(w)} = -2 



(1) 



is satisfied, where ep and em are the dielectric functions 
of the metal particle and molecular environment. Indeed, 
even in this limit, Eq. ^ is difficult to solve when both 
the metal and environment exhibit resonant behavior in 
the same spectral regime. 

We calibrate our first numerical description of LSPR 
spectroscopy through an approximate modeling of a 
series of recent Van Duyne experiments [3, EBl prob- 
ing the influence of electronically-resonant molecules 
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upon the LSPR of bare silver nanotriangles. Em- 
phasis is placed not upon the exact modeling of any 
one of these experiments but rather upon under- 
standing behavior generic to all. To this end, we 
choose to replace the reported [2,3,7,8,12,13,17,18- 
octakis(propyl)porphyrazinato]magnesium(II), iron(II) 
tris-2,2'-bipyridine, and rhodamine 6G with the simpler 
pyridine molecule, treated fully quantum mechanically, 
and to approximate the silver nanotriangles by classical 
silver nanospheroids, which also have tunable aspect 
ratios and resonances; as in the experiments, this 
tunability will allow us to tune the LSPR into or detune 
away from one or more of the molecular resonances. 
To our knowledge there are no previous first-principles 
models of plasmon-enhanced spectroscopy that have 
investigated the self-consistent coupling between molec- 
ular and true plasmonic resonances within the context 
of LSPR spectroscopy for real systems like this. And it 
is a main point of this paper to report the achievement 
of such a computation. 

By first principles we mean a theoretical description 
of the molecule-metal system in which the mutual cou- 
pling between molecule and metal and its effect on 
the molecular-electronic structure and electromagnetic- 
scattering properties of the combined target is computed 
without varying parameters to fit experiment, while the 
underlying response of the metal is completely governed 
by the phenomenological dielectric data for bulk metals; 
see, e.g., Ref. jjJi]. Other noteworthy theoretical ap- 
proaches (some of which are first principles) have been 
applied to a variety of plasmon-enhanced processes such 
surface-enhanced Raman [J, 0, Q and fluorescence spec- 
troscopy [3| as well as the transfer of energy through ar- 
rays of metal nanoscale particles Q. For example, Corni 
and Tomasi extend the polarizable-continuum model of 
molecular solvation to treat the interaction of a gen- 
eral molecular-electronic system with a nearby classi- 
cal polarizable metal particle and compute Raman en- 
hancement factors as a function of excitation energy [l[ . 
While Jensen et al. compute SERS spectra from a va- 
riety of molecules adsorbed onto small metal clusters, 
intended to mimic nanoscale metal particles, with mod- 
ern computational-chemistry methods, but are forced to 
empirically parametrize the polarization interaction be- 
tween the molecule and the remainder of the metal Q. 
Johansson et al. describe both SERS and fluorescence 
of a two-state molecular-electronic system placed in the 
junction between two metal spheres within a density- 
matrix formalism that phenomenologically includes both 
radiative and nonradiative relaxation mechanisms 0. 
And, within the context of energy transfer, Neuhauser 
and Lopata combine a finite-difference time-domain elec- 
trodynamics description of the metal particle and ex- 
ternal field with a random-phase approximation of a 
two-state molecular-electronic system's explicitly time- 
dependent dynamics, together with a phenomenological 



damping of the molecular excitations due to the pres- 
ence of the metal Lastly, Masiello and Schatz rig- 
orously combine modern molecular electronic-structure 
theory with a continuum-electrodynamics description of 
the metal, including the self-consistent coupling between 
the two, to enable practical computation of SERS in real 
systems From a more general perspective, we also 
note the recent work of Yuen-Zhou et al. to extend the 
Runge-Gross theorem to open quantum systems, thereby 
treating the interactions between molecular electrons as 
well as their coupling to a class of Markovian and non- 
Markovian environments on equal footing Q . 

Within our formalism, care is taken to self-consistently 
incorporate the polarization effects between a general 
quantum-mechanical molecular system and an arbitrary 
nanoscale metal structure supporting classical plasmon 
modes of excitation in the presence of the classical elec- 
tromagnetic field. Generalizing and extending the orig- 
inal work of Gersten and Nitzan [l^, the dynamics of 
the molecular adsorbate is built upon modern electronic- 
structure quantum-chemical theory and is simultaneously 
coupled to a full computational-electrodynamics descrip- 
tion of the metallic plasmons and their associated elec- 
tromagnetic fields. Symbolically, the linear response of 
such an interacting many-electron molecule to an exter- 
nal electromagnetic potential is described by the in- 
duced electronic density 

p(i)(x,i) = (l/ri)/d3.x'dt'n«(x,i;x',i')(-e)$i;(x',t'), 

(2) 

where the plasmonic polarization effects are encapsulated 
within the molecular polarization propagator U. Solution 
of this microscopic equation forms a large part of the 
following discussion, as does its self-consistent incorpo- 
ration within a macroscopic-electrodynamics description 
of the combined molecule-metal system's electromagnetic 
scattering. A related, but purely theoretical development 
was previously reported within the context of surface- 
enhanced Raman scattering Q; however, there we con- 
sidered the inelastic scattering of light from a molecular 
target, while here we consider the elastic scattering of 
light from nanoscale metal particles as well as the re- 
sponse of the surrounding molecular environment, in- 
cluding their mutual coupling. Yet since much of the 
formalism is common, here we choose to place emphasis 
upon the numerical implementation and refer the inter- 
ested reader to Ref. [4|] for more theoretical details. As in 
our previous work [J| , chemical- interaction effects where 
the molecular wave function would extend onto the metal 
and vice versa are neglected. 

Basic to our development is the dequantization of both 
the electromagnetic field (E,B) (E,B) and the plas- 
mon field fl — *■ ^/n in their large occupation number 
limits, leaving only the molecular system quantum me- 
chanical. A Born-Oppenheimer separation of electronic 
and nuclear degrees of freedom in the molecule facilitates 
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the use of many-body perturbation theory to couple in 
the electronic polarization effects of the metal and elec- 
tromagnetic field. Acting as a second polarizable body 
with plasmonic density n, the metal, if allowed to in- 
teract with the molecular system with electronic density 
p through the potential Vp{t) — —ej d'^xp{yi,t)^p{x,t) 
with $p(x, i) — —e J d'^a:'n(x', i)/|x — x'|, will generate 
a proper self-energy correction S* to the noninteracting 
one-body Green's function G of the isolated molecule. 
Truncation of E* « E*^ at second order in the den- 
sity Sn = n — Uq induced in the metal by the molecule 
allows for a partial resummation of this polarization in- 
teraction through the inversion of the one-body Green's 
function = — S*^ of the interacting molecule. 

It is through the computation of {G^^)^^ that the two 



bodies repeatedly polarize each other, to all orders in 
perturbation theory, thereby allowing a self consistency 
to be reached between the molecule and its image as seen 
in the metal. 



We base this formalism upon an initial mean-field de- 
scription of the isolated molecular system at either the 
Hartree-Fock (HF) or Kohn-Sham density-functional the- 
ory (DFT) level, and subsequently construct the matrix 
elements of and in that basis. Adopting an 
isotropic Lorentz model of the metal's complex-valued 
polarizability cxp allows for the analytic evaluation of 
the convolution appearing within the second-order time- 
ordered self energy 



duj' 



h[E*^]pq{uj) = ^{p\ - ex|r) • A • / —ap{uj')Grs{Lo + lu')-A - {s\ - ex|q) 
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(3) 



with ujq — Hq — (70/2)^. It is expressed in terms of 
the noninteracting molecular orbitals {</>}, orbital ener- 
gies {e*^}, and occupation numbers {p'^} of the nonin- 
teracting molecule at its ground-state equilibrium nu- 
clear geometry; these quantities are readily obtained 
from any standard electronic-structure computation. As 
usual, indices i,j,k,l,... {a,b,c,d,...) refer to orbitals 
that are occupied (unoccupied) in the reference state 
(either mean- field or Dyson), while indices p, q,r, s, . . . 
are arbitrary. Unless otherwise noted, all integrals are 
taken over the full range of coordinates and all sums run 
over the complete model space. And, as in Ref. [3], 
A^"" = {SR^R" - S^'')/R^ with X,a = x,y,z are matrix 
elements of the instantaneous near-zone dipole tensor A 
with R = i?R the vector connecting molecule and image 
[20| : see Fig. [T] Since it is not our intent to compute 
properties intrinsic to the metal, the parameters flo, 70, 
and ao representing its resonance positions, line widths, 
and amplitudes are chosen to fit the polarizabilities of 
nanoscale metallic structures within some model, e.g., 
Clausius-Mossotti or Gans for spheres or spheroids, or 
the lattice-dispersion relation [2l| for arbitrarily shaped 
targets, and depend upon its underlying dielectric re- 
sponse Extension to multiple Lorentz oscillators of 
the form ap{uj) ~ ap(oo) — 7 + ii^Jj — ^^7) (or 

even multiple Drude and/or free-electron gas oscillators) 



is straightforward as '^(^2) linear in ap. 

Under the assumption that the self energy varies 
slowly with frequency, diagonalization of the complex- 
symmetric matrix ^"^(/lo;) « G^^(ftw)— S*^(e°) at each 
frequency lu 

Qo^ihw) = v--^{huj)g-^{nuj)\]{nw) = hu}^s (4) 

results in a set of interacting and complex-valued Dyson 
orbitals {x} making up, e.g., the columns of the right 
eigenvector matrix U and corresponding orbital energies 
{e} = {e° + /ip* ]75(e°)} - {£« + fiA(eO) + ihTis'')} 
forming the matrix elements of the diagonal matrix 
hij — e. From Eq. (jH), the diagonal and interacting one- 
body molecular-electronic Green's function has the time- 
ordered form 

[SoXqif^) = ^P9[(i-Pp)(^-£p+«4)+/'p(^-4-^4)]' 

(5) 

where both £'' — e'^ + ftA(£°) and = fir(e°) are ap- 
proximated as frequency independent, but do vary from 
single-particle state to single-particle state. It is the 
imaginary part of e that accounts for the finite lifetime 
h/e^ of the electronically excited states of the molecule 
due to their damping by the polarization (image) inter- 
action with the metal, while the real part is responsible 
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for the corresponding shift of molecular energies. Con- 
struction and diagonalization of the full self-energy ma- 
trix SS*^ allows both of the matrices HA. and hT to be 
extracted. 

Equipped with these Dyson orbitals and orbital en- 
ergies, the density response Sp — p — (p) of the in- 
teracting molecule to an external perturbation VE(t) = 
—e J d^xp{x, t)<i>£;(x, t) can be computed from the causal 
response kernel 11^ known as the retarded polarization 
propagator. The lowest order (linear) response of the 
density is (Sp{y:,t))E = p''^\^,t)-\ , where the expec- 
tation (•••)£; is taken in the interacting state (• • • ) ad- 
ditionally perturbed by Ve- It has the matrix elements 
piWuj) = (l/ft)E..n«,,(c^)(-e)$f,H. Once n« and 



Ve are known, then all observable properties related to 
the linear response of the system can be computed. For 
example, the linear induced-dipole moment d^^^{t) = 
J dt'aM{t, t') ■ E(t') is obtained from Eq. ^ [or its ma- 
trix elements p^pq] by multiplication by J (fix{—ex.) [or by 
X]pq(p| "'S^l'i')] the electric-dipole interaction approxi- 
mation with <i>£;(x, t) = — x-E(t). All that remains is the 
computation of 11, whose matrix elements are related to 
the interacting one-body Green's function Q through the 
convolution npgrs(w) — J{dLu' /2Tri)Gpri^ + ^')Gsqi^')- 
We analytically perform this integral by first rotating Q 
into its diagonal frame where it takes on the simple time- 
ordered form (O and then extract its retarded component 



T^pqrsi^)^ E Upp-U,s'[Tl%,q'r's'ii0)U-lU-l 
p' q'r' s' 



- Uvp-Uss-5p.r'5s'q- y^l'^.yj^^^^q.q-^r'r 

p'q'r's' y P ^ H P 



(6) 



under the assumption that Upq{uj) ~ Upq{eq) = Upq. Rotation back to the original frame is subsequently performed 
numerically. By multiplication of in^(x, t; x', t') = <d{t — t'){[p{x,t), p{x' ,t')]) with / d^xd'^x' {—e'x.){—ex') or con- 
traction of p^J with (p| — ex|g), the one-body polarizability a]\{ of the molecule including its polarization (image) 
interaction with a nearby plasmon supporting metal is found to be 



^Y.(^\-ex^\p)U^^^,{uj){r\-ex^\s) 



{a\d^\i){i\d''\a) 



{a\d''\i){i\d^\a) 



(£- - £[) + tiel + el) + hLu (£- - <) - i{el + <) 
I 



(7) 



with one-body dipole-matrix elements (pld^jf?) = 
'I2rs^ps^i^\ ~ ex^\r)Urq rotated into the Dyson basis 
denoted by the symbol | ). This expression generalizes 
the standard one-body approximation to the Kramers- 
Heisenberg polarizability and further extends the work 
of Norman et al. [2^ and Jensen et al. [1^, where 
the infinitesimals 0^ appearing in the denominators of 
the molecular polarizability are replaced by an empiri- 
cal parameter SFoxp that mimics the molecule-metal in- 
teraction and defines the polarizability on resonance, to 
explicitly and rigorously include the resonance widths 
£^ + e]. It is, of course, important to note that Eq. ^ 
reduces to the usual one-body molecular polarizability 
of an isolated molecule in the limit where 
£* = 0+, i.e., when Vp — > 0. In that case U 
ip\d^\q) {p\-ex^\q). 



-0 



and 
1 and 



To go beyond the mean-field approximation of the 
molecular electron-electron interactions in the calcu- 
lation of the ckM, perturbation theory may be ad- 



ditionally performed on the time-ordered polarization 
propagator 11. Namely, the two-body Coulomb po- 
tential VAiit) = J d^xd^x'p{K,t)VM{y^,x')p{x',t) ^ 
(e^/2) / d^xd^x'p{x,t)p{x.',t)/\x — x'\ generates a proper 
self-energy correction K* to the interacting 11 of the 
molecule-metal system. Truncation of K* « K 



(1) 



at 



first order in perturbation theory allows for a partial re- 
summation of the effects of electron-electron interaction 
at the level of the random-phase approximation (RPA), 
through the inversion of the RPA polarization propaga- 



tor IIt; 



n- 



^RPA — ~ '^(l)' ^'^ 

the poles of IIrpa or the zero eigenvalues of IIj^p^ that 
the low-amplitude collective excitations of the molecular- 
electronic system are located. These zeroes can be found 
by solving the associated pseudoeigenvalue problem 
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the Dyson 

■'q'r's' •~'pp"~'ss' [n 



basis, where [IIrpa ] 

RPAJp'i3'r's"-^g'q ^ ' 
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and 
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n[K* 



(i)ipqr 



-ips\VM\rq) + {sp\VM\rq). Eqs. 



may be decomposed into {i, a) and (a, i) components 



+ [{ba\VM\3i) - {ab\VMmp^^^{u^)] = 

3b 



(9) 



r 



where they generahze the well-known RPA equations 
to complex-valued orbitals and orbital energies. In the 
noninteracting limit, i.e., when Vp — > 0, the standard 
RPA equations are recovered with Aiajb — <5ij(5fca(£° — 
e^) -I- {ia\VM\bi) - {aj\VM\bi) and Biajb = {ha\VM\ii) - 
(a6|VA/|jj). Solution of these generalized RPA equations 
([9|) followed by a rotation back to the original basis 
results in new eigenenergies that account for all pos- 
sible single-electron excitations and single-electron de- 
excitations from the interacting molecular reference sys- 
tem in the Condon approximation; setting _B = would 
be equivalent to performing a configuration interaction 
with all possible single excitations from the mean-field 
reference, also called the Tamm-Dancoff approximation. 
The resulting polarizability —ha.^^j^{Lu) = X]pqrs(9l ~ 
ex|p)([n^p^]~^)pg^^j,((i;)(r| — ex|s) built from these ener- 
gies would improve upon that in Eq. ^ to include this 
limited class of electronic excitations. 

We numerically implement these molecular interac- 
tions with the metal and the external field due to the per- 
turbations Vp and Ve as a post-processing step within a 
modified version of the Q-Chem quantum-chemistry soft- 
ware package [i^, but do not treat any RPA- level elec- 
tronic excitations induced by Vm at this time. These for- 
mer interactions are built upon an initial mean-field cal- 
culation of the molecule's electronic ground state, which, 
for the examples presented in the following, will be DPT. 
Dyson orbitals and corresponding orbital energies are not 
brought to self consistency in the sense that the self en- 
ergy ([3]) and interacting one-body Green's function dH) 
are computed directly and only once from the nonin- 
teracting Kohn-Sham orbitals and orbital energies {0} 
and {e"}. Continuum-electromagnetic scattering prop- 
erties of the combined metal-molecule system are sub- 
sequently computed through a version of the discrete- 
dipole approximation (DDA) [2lj that has been modified 
to directly accept the microscopic and isotropic molecu- 
lar polarizability as input rather than the molecular sys- 
tem's effective macroscopic index of refraction or dielec- 
tric function. 

As a first application of our formalism, we compute 



the linear response and scattering of pyridine interact- 
ing with a nearby nanoscale silver spheroidal particle; 
see Fig. [2j By deforming the spheroid's aspect ratio 
from prolate (left panel) with radii (50,100,50) nm, to 
sphere (middle panel) with radii (50,50,50) nm, to oblate 
(right panel) with radii (50,1,50) nm, the particle's y- 
axis plasmon resonance may be tuned across a portion 
of the electronic spectrum of the pyridine-nanospheroid 
system, thereby allowing us to systematically explore 
their interaction. Sensitivity of the LSPR position to 
the surrounding environment's refractive index in the 
monolayer-coverage limit of a variety of electronically 
resonant molecules on silver nanotriangles has recently 
been observed 



14 



15l | and serves to motivate our first- 
principles calculations. 

Our computations involve two steps: first, the micro- 
scopic quantum-mechanical response of a single pyridine 
molecule is computed as it interacts with its image in 
a nearby silver nanoscale spheroidal particle, described 
analytically by either Clausius-Mossotti or Cans mod- 
els for ap(ep); see Fig. 2, upper panels. Second, the 
resulting isotropic pyridine polarizability is incorporated 
into a DDA description of the macroscopic electromag- 
netic scattering of a monolayer-covered silver spheroid, 
where the monolayer and metal are discretized into a 
collection of polarizable points with polarizabilities de- 
fined from step one for the molecule and by the lattice- 
dispersion relation ,21i] for the metal; see Fig. 2, lower 
panels. In both steps the underlying dielectric data for 
silver is taken from Ref. [l^. 

To further elaborate, first, we fit multiple Lorcntz oscil- 
lators to CKp, which together with a DFT-VWN ground- 
state electronic-structure calculation in a 6-31g* basis, 
define the self energy S*^ ([3]), one-body Green's function 
Q and polarization propagator II-'^ and polarizability 
a A/ of the interacting molecule. From the imaginary 
component of olmi we display the averaged extinction 
cross section of a single pyridine molecule as a func- 
tion of energy and spheroid aspect ratio in Fig. [2] (upper, 
blue), together with that computed (upper, black) using 
the approximations of Refs. 2^, |23|; the latter chosen 
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FIG. 2: Averaged extinction cross sections cr*xt = 4iTklm[a%f + + Q:|f]/3 (upper) and extinction coefficients (Tcxt/(cross- 
sectional area) (lower) of the coupled pyridine-silver nanospheroid system as a function of excitation energy and plasmon- 
resonance position. Explicit treatment of the repeated polarization interaction between a single pyridine molecule and its 
image in a nearby metal spheroidal particle reveals that the pyridine polarizability and extinction (upper, blue) are sensitive 
to nanoparticle shape, LSPR position and width, and relative molecule-metal separation and orientation, here chosen so that 
molecule and image are coUinear with each other {(f = 0) and with the field polarization direction, and separated by i? = 1.6 nm. 
For comparison, we also compute the pyridine extinction cross section based upon an empirical and constant parametrization 
of the molecule-metal interaction E* ~ iFexp (upper, black) that is consistent with the work in Ref. Progression from a 

pyridine-coated oblate spheroid (right) with radii (50,1,50) nm, to sphere (middle) with radii (50,50,50) nm, to prolate spheroid 
(left) with radii (50,100,50) nm allows for the y-axis plasmon resonance to be tuned from nearly resonant with to red-detuned 
from the lowest-lying molecular resonance of the interacting molecule. Extinction coefficients of the coupled system in the 
monolayer-coverage limit (lower, red) and bare metal (lower, black) are shown with the excitation field polarized along the y 
axis as well as polarized along the z axis in the oblate case for the bare metal only (lower right, dotted black). Insets not drawn 
to scale. 



to empirically parametrize the molecule-metal interac- 
tion S* w iFoxp directly by a line width of ^iFcxp = 0.1 
eV that is independent of molecular state, energy, rela- 
tive molecule-metal orientation and separation, and plas- 
monic structure. To compare both approaches, we orient 
pyridine's dipole moment to be normal to the metallic 
surface (a choice of = in the representation shown in 
Fig. [T|), along the field-polarization axis, and separated 
from its image by i? = 1.6 nm; this distance is chosen be- 
cause it falls in the middle of a regime, from ~ 1 — 3 nm, 
where image effects seem to be important, i.e., are com- 
parable in magnitude to the homogeneous broadening, 
while smaller separations reveal a (expected) breakdown 
of our perturbative approach. We also adopt a reason- 
able and constant value for the homogeneous line width of 
pyridine to be hjH = 0.05 eV based on analogy with es- 
timates for other molecules [iili^,!!!]. Since our model 
only describes the modification of the molecule's dynam- 
ics due to its interaction with the plasmon and does not 
include any intrinsic- or homogeneous-broadening mech- 



anisms it is necessary to add H'jh before comparing with 
experiment; however, we assume that such mechanisms 
are completely decoupled from the relevant plasmon en- 
hancement, thereby justifying such a separation. In Fig. 
[21 it is important to note that the shifting and broad- 
ening of molecular-excited states due to their interac- 
tion with the metal are dependent not only upon the 
LSPR position of the metal but also upon its strength; 
e.g., the lowest- lying pyridine excited state is nearly res- 
onant with the oblate spheroid's y-axis dipole-plasmon 
resonance, yet its interaction with the oblate spheroid is 
not as strong as with the red-detuned sphere or prolate 
spheroid because of the negligible strength of the oblate 
spheroid's LSPR; note the change in magnitude of the 
extinction coefficient at the LSPR maximum of the pro- 
late spheroid, sphere, and oblate spheroid in the lower 
three panels of Fig. (5] Maximizing the two, oscillator 
strength and spectral overlap, would lead to the greatest 
interaction effects modifying the molecular polarizability. 

Second, we add a pyridine monolayer to a DDA-level 
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continuum-electrodynamics calculation of the bare sil- 
ver spheroid, where the monolayer and metal are built 
from a collection of polarizable points with isotropic po- 
larizabilities ctM computed from step one and ap(ep) 
computed from the lattice dispersion relation [2lj; each 
DDA point occupies a cubic volume of (2 nm)^ with 
other choices of nearby dipole spacing yielding similar re- 
sults. Here we point out that while a.M accounts for the 
repeated polarization interaction between molecule and 
metal, it does not include any direct quantum-mechanical 
molecule-molecule couplings; however, it would be inter- 
esting to explore the effect of different molecular aggre- 
gation schemes [3] and how they interact with the metal 
with our methodology in the future. Numerical solution 
of this combined DDA accounts for the self-consistent 
electric-dipole interaction between all of the metallic and 
molecular dipoles among themselves and, additionally, 
with each other. The resulting extinction coefficients 
Cext/ (cross-sectional area) of the combined and coupled 
system (lower, red) as well as those of the bare silver 
spheroid (lower, black) are displayed in Fig. [2] as a func- 
tion of excitation energy and spheroid aspect ratio. The 
progression from no shift (AA — nm, oblate) to a red 
shift (sphere, AA — 7.88 nm; prolate, AA = 26.02 nm) 
of the combined system's LSPR wavelength A with re- 
spect to that of the bare metal, as the plasmonic and 
lowest-lying molecular resonances are detuned, is consis- 
tent with the observations reported in Refs. [lJ,[lB|- Also 
note that while the pyridine monolayer does not affect the 
position of the oblate spheroid's y-axis dipole-plasmon 
resonance, it does mediate the coupling of incident light 
to the previously dark z-axis plasmon resonance; in par- 
ticular, the feature appearing at 2.6 eV (lower right, red) 
should be compared to the new plasmon resonance at 2.7 
eV of the bare oblate spheroid excited with z-polarized 
light (lower right, dotted black). 

As a further test of our formalism we examine the fic- 
titious interaction of pyridine with a 20-nm radius quasi- 
silver sphere, whose underlying dielectric function [18| is 
shifted, by hand, to allow its dipole-plasmon resonance to 
be tunable across several low-lying molecular-electronic 
excitations. (This is normally not possible for pyridine- 
silver nanoparticle system.) By varying this dielectric 
shift (DS), each bullet in Fig. [3] (connected in red) repre- 
sents a separate computation of the extinction spectrum 
and resulting plasmon-resonance shift (PRS) of the com- 
bined molecule-quasi-metal system with respect to that 
of the bare quasi-metal; nm DS corresponds to the real 
pyridine-silver system with an LSPR wavelength of 370 
nm. For comparison, the PRS of a 20-nm radius quasi- 
silver nanosphere coated with a uniform and nonresonant 
layer with index of refraction n = 1.5 is shown as the 
dashed black curve. Superimposed on the plot is the ex- 
tinction cross section of a single pyridine molecule (blue) 
interacting with its image, located i? = 1.6 nm away and 
with dipole moment aligned with the field-polarization 
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FIG. 3: Plasmon-resonance shift of a pyridine-coated 20- 
nm radius quasi-silver nanosphere as a function of the bare 
sphere's tunable plasmon-resonance position (black bullets, 
connected in red). Variation, in increments of 10 nm, of the 
sphere's underlying dielectric response (from -130 nm DS to 
+60 nm DS) with respect to that of silver (0 nm DS taken 
from Ref. J8]) allows us to tune the bare quasi-metal's LSPR 
across several low-lying molecular-electronic excitations. For 
comparison the PRS of a 20-nm radius quasi-silver nanosphere 
coated with a nonresonant monolayer with index of refraction 
n = 1.5 is displayed as the dashed black curve. Correspond- 
ing to a nm DS, the computed extinction spectrum of a sin- 
gle pyridine molecule interacting with its image in the silver 
nanosphere (blue) is superimposed for reference along with 
that computed from an empirical treatment of the pyridine- 
metal interaction (black) consistent with Ref. [23| . 



direction (a choice of 1^9 = in the representation shown 
in Fig. [1]), in the silver sphere (for the case having 
nm DS), as well as that computed from an empirical 
parametrization of the molecule-metal interaction (black) 
consistent with Ref. [2^. (These curves are analogous 
to the blue and black curves in the upper panels in Fig. 
[21) Our computations reveal that the PRS increases as 
the DS becomes more negative. This is the expected re- 
sult for a nonresonant scatterer interacting with a metal 
particle, as the shift of the plasmon resonance should 
become larger as the particle becomes more polarizable, 
and in the present application this corresponds to going 
to the blue. However, we also notice dips in the PRS 
at ~315 nm and ~245 nm, which coincide with molec- 
ular resonances in pyridine (blue curve). The dips are 
understood to arise because the localized but negative 
(anomalous) polarizability of pyridine destructively in- 
terferes with the positive (normal) portion of that of the 
particle as the LSPR passes through resonance [13, [3] ■ 
This behavior is related to behavior found in the exper- 
iments in Refs. 



14l . Il5j , which show a pronounced dip 



in the PRS as the plasmon is tuned through the molec- 
ular resonance. In this case, the dips are smaller due to 
the weakness of the pyridine transitions, but the physical 
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effect is the same. 

In conclusion, we present the first numerical implemen- 
tation of our general many-body description of a variety 
of plasmon-enhanced linear spectroscopies. We briefly re- 
view our theoretical formalism in order to motivate and 
detail our computational efforts in linking together both 
Q-Chem [2^ and the DDA 21 1 to self-consistently and 
practically describe the linear-response and scattering 
properties of an arbitrary quantum-mechanical molecule 
interacting with a nearby classical nanoscale metal in the 
presence of an external continuum electromagnetic field. 
As a preliminary test of our methodology we compute 
the linear response and scattering of the pyridine-silver 
nanospheroid system as a function of excitation energy 
and spheroid aspect ratio. Variation of the latter allows 
us to tunc the nanoparticle's LSPR position from nearly 
resonant with to red-detuned from the molecule's lowest- 
lying electronic excitation; the positions and widths of 
such are here computed from first principles. Interac- 
tions between these resonances are highlighted and qual- 
itatively discussed in connection to recent experimental 
observations of similar systems 14 , 1^ . In addition, we 



examine the fictitious interaction between pyridine and 
a quasi-silver nanosphere with tunable LSPR position in 
order to further explore the coupling of the plasmon to 
multiple electronic excitations in the molecule. Quan- 
titative and rigorous modeling of a variety of plasmon- 
enhanced linear elastic and inelastic spectroscopies is now 
enabled with our approach and will serve to motivate our 
future investigations. 
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